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ABSTRACT 

A study is made of the behavior of massive black holes in disk galaxies that have received an 
impulsive kick from a merger or a sustained acceleration from an asymmetric jet. The motion 
of the gas, stars, dark matter, and massive black hole are calculated using the GADGET-2 
simulation code. The massive black hole escapes the galaxy for kick velocities above about 
600 km s -1 or accelerations above about 4 x 10~ 8 cm s~ 2 over time-scales of the order of 10 s yr. 
For smaller velocity kicks or smaller accelerations, the black hole oscillates about the center of 
mass with a frequency which decreases as the kick velocity or acceleration increases. The black 
hole displacements may give rise to observable nonaxisymmetries in the morphology and dynamics 
of the stellar and gaseous disk of the galaxy. In some cases the dynamical center of the galaxy is 
seen to be displaced towards the direction of the BH acceleration with a characteristic "tongue-" 
shaped extension of the velocity contours on the side of the galaxy opposite the acceleration. 

Subject headings: galaxies: kinematics and dynamics — galaxies: structure — galaxies: nuclei 



1. Introduction 

Recent breakthroughs in numerical General 
Relativity have led to predictions of large recoil ve- 
locities of merged binary black-holes (BHs) as the 
binary radiates away linear momentum as grav- 
i tational waves during the final stages of merge r 



the merging BHs and a large ratio of the mass es 



of the two initial BHs ( Bogdanovic et al. 2007t h 



In another situation, the asymmetric or one- 
sided jet from the disk around a massive BH can 
give a steady acceleration of th e black hole o ver th e 



([Gonzalez et al. I (|2007l ) ICampanelli et al. I(|2007l )) 



Salpeter time-s cale ~ 10 8 yr ([Shklovskv I (|1982| ); 



Tsvganl(l2007l)) and this may p ush the BH out of 



The velocity or "kick" of the merged black-hole 
in the nucleus of a galaxy can be as large as 
4000 km s , which would be more than enough 
to eject the BH from the galaxy. However, it is 
more likely that typical kick velocities are much 
smaller (< 200 km s^ 1 ) owing to gas accretion by 



the galaxy. Wang et al. I (j 1992h pointed out that 
some combinations of dipole and quadrupole sym- 
metry magnetic field components threading the 
accretion disk of the BH can give a much stronger 
magnetically driven jet on one side of the disk than 
on the other. It is of interest that a bright quasar 
has been discovered which is not inside a massive 
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galaxy (jMagain et al. 1120051 ) and so it may have 



been ejected from a galaxy. 

Field galaxies themselves may be morpholog- 
ically or dynamically nonaxisym metric even in 
the absence of tidal companions (jBaldwin et a 
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1980). Up to 30% of field galaxi e s are nonaxisym- 



metric (jKornreich et al. I (|1998l ). IZaritskv fc Rix 
(|1997l )). If central supermassive BHs receive im- 
pulses from mergers or asymmetric jets, causing 
motion in the host galaxy, some of that motion 
may be transmitted to the galaxy via dynamical 
friction. We may expect that in such galaxies, 
morphological or dynamical effects may be seen in 
the stellar or gas distribution as a result of this 
transfer of momentum. 

The objective of this work is to determine the 
BH displacement as a function of time for different 
initial kick velocities and directions, different ac- 
celerations due to a one sided jet, and different BH 
masses. A second objective is to study influence 
the BH displacement has on the morphology and 
the kinematics of the host galaxy and determine 
whether nonaxisymmetries in disk galaxies are at- 
tributable to motions of the BH. We have chosen 
to concentrate on the morphology of disk galax- 
ies rather than ellipticals because the dynamics 
of disk galaxies, being supported by bulk rotation 
rather than random motions, are more likely to 
exhibit signs of disturbance if the motion of the 
BH interferes with that rotation. 

Section [2] of the paper describes the initial con- 
ditions and the numerical methods used. Section[3] 
discusses the results, first for the case of velocity 
kicks of a 10 8 M Q BH, and then of a 1O 9 M BH. 
Section 0] discusses the results for cases of one- 
sided jets, first for the 1O 8 M BH, and then for 
the 1O 9 M BH. Section [5] gives the conclusions of 
this work. 

2. Numerical Methods and Initial Condi- 
tions 

The computations for this experiment were per- 
formed using the N-body, Smoothed Pa rticle Hy- 



drody namics (SPH) GADGET-2 code of lSpringel 
(|2005l ) . We have added the ability to assign a par- 



ticular particle an arbitrary acceleration, in addi- 
tion to the calculated gravitational acceleration, in 
order to simulate forces due to an asymmetric jet. 
This particle is given significant mass and placed 
initially at the center of the galaxy. 

The considered galaxy consists of a dark halo, 
a stellar disk and bulge, and a gas disk. The 
gas particles are treated gravitationally and hy- 
drodynamically, while the other particles interact 



gravitationally. The halo, disk, and bulge parti- 
cles differ in mass and spatial distribution, but are 
otherwise treated identically in the co mputation. 



The ce ntral stellar bulge is modeled as a lPlummer 



(|1915l ) sphere, which has a density distribution of 
a relaxed n = 5 polytrope given by: 



p(r) 
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(1) 



with bulge scale radius Rb = 1 kpc and mass 
M B = 1O 1O M . 

This bulge is placed at the center of a softened 
isothermal dark matter halo of density profile 



p(r) 



1 + (r/r c f 



(2) 



total mass 5.6 x lO n M0, with a cutoff radius 
50 kpc, and core radius r c — 2 kpc. A stellar and 
gas disk of scale radius tq = 3.5 kpc and cutoff 
radius R = 50 kpc of exponential surface density 



E(r) 



M D 



2-ktI [1 - e~ R / r ° (l+R/r )] 
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(3) 



surrounds the bulge with total mass Md = 2.8 x 
1O 1O M0, at a gas mass fraction of 25%. In the ver- 
tical direction, the disk particles are given velocity 
dispersions and uniformly random initial positions 
corresponding to a disk of 1 kpc thickness. 

The total dynamical mass of the galaxy model 
is therefore M tot = 5.98 x lO n M . This model 
was chosen to conform to parameters observed 
in giant disk galaxie s su ch as those described by 
Rownd et arT(|l994 ) and lDickev et al. I (|l990h . in 



particular, a circular velocity of 220 km/s with a 
10 km/s velocity dispersion in the disk. The total 
number of simulation particles used in the model 
was 2 x 10 6 for most model runs. 

In one sequence of runs, we gave the central 
black hole an initial impulsive velocity "kick" to 
simulate the recoil from the merging of two black 
holes. In one sequence, the kick was in the plane 
of the galaxy, and in the other, it was normal to 
the galaxy plane. The black hole was given no 
further acceleration. In another sequence of runs, 
the BH was given a steady acceleration of order 
10 -8 cm s _1 for a time of the order of the Salpeter 
time in order to simulate the influence of a one- 
sided jet. Each of these sequences was run for 
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a BH mass of 10 8 M Q and a larger black hole of 
10 9 A/ Q . 

During each run, the position, velocity, and ac- 
celeration of the BH was output for each timestep. 
Data dumps containing positions, velocities, ac- 
celerations, and (for the gas particles) tempera- 
tures were generated periodically, usually every 
0.01 Gyr of simulation time. These data were read 
into the IDL data processing environment, where 
it was possible to display the results as seen from a 
variety of viewing angles. For a given viewing an- 
gle, gas particles were convolved with a Gaussian 
simulated telescope beam and displayed as loga- 
rithmic intensity contours overlaid on a diagram 
of stellar particle positions. Gas velocities were 
similarly convolved and used to produce radial ve- 
locity maps. 

Additional sample runs were conducted with 
isothermal halo core radius 0.5 kpc, with no 



radius, and with the 



of 



Navarro et al. 



halo density profile 
In each of these cases 



( 19971) . 

the results were similar and the details of the be- 
havior did not substantially depend on the nature 
of the halo density profile. 

3. Velocity "Kicks" 

3.1. Black Hole Mass 1O 8 M 

In this series of model runs, the BH was placed 
near the center of mass of the galaxy and given 
an initial velocity in the plane of the disk. The 
series of initial velocities started at zero and pro- 
gressed in increments of 50 km s" 1 up to a velocity 
of 650 km s _1 , for which the central BH escaped 
the galaxy. It is interesting to note that this crit- 
ical kick velocity is significantly larger than the 
escape velocity of the center of the galaxy, which 
was measured to be 450 km s _1 . We attribute this 
excess to dynamical friction. 

The general behavior of the BH after receiv- 
ing the kick was to oscillate about the model cen- 
ter of mass approximately sinusoidally (Figure [I]). 
The amplitude and frequency of the motion re- 
main roughly constant over a period of 0.5 Gyr, 
but depend on the kick magnitude as shown in 
Figure [21 Error bars represent the la standard 
deviation of these values as calculated for each cy- 
cle. For kick velocities below 400 km s _1 , the 
relation between amplitude and initial velocity 
V is linear, with A(V )/kpc = (9.30 ± 0.53) x 



lO^Vb/kms -1 + (0.088 ± 0.082). For initial ve- 
locities above 400 km s _1 , the BH amplitude is 
beyond the visible disk and becomes exponentially 
large until it escapes at 650 km s -1 . The frequency 
of the oscillation is also linear with /(Vb)/Gyr _1 = 
(17.64 ± 0.57) - (0.0317 ± 0.0015) F /km s -1 . 

Smaller kicks resulted in the BH remaining in 
the more massive parts of the galaxy within 10 kpc 
for their entire cycles, allowing momentum from 
the kick to be efficiently transferred to the galaxy 
via dynamical friction. In these cases, damping 
of the oscillation was on timescales of 0.7 Gyr. 
The galaxy center of mass thus gained velocity in 
the direction of kick, resulting in the offsets de- 
picted in the low kick amplitude curves in Fig- 
ure [2] Larger kicks, while imparting greater mo- 
mentum to the BH, caused it to spend a great deal 
of time in the thin regions of the galaxy beyond 
r = 10 kpc. This greatly reduces the efficiency 
of the dynamical friction, lengthening damping 
timescales to 2.0 Gyr and reducing the momen- 
tum transferred to the galaxy. 

With the frequency of the BH motion to = 2nf 
known as a function of the amplitude of its motion, 
we can derive the effective potential for its motion, 
L/ e ff. For a one-dimensional oscillator, 



R = 



dE 



V2 Jo uj(E)^U cff - E 



(4) 



( Landau fc Lifshitz I ll960h . A first approximation 
to our simulation results gives lu — luq(1 — Vo/V m ), 
where Vq is the BH kick velocity such that E = 
V 2 /2, and where uiq and V m are constants charac- 
teristic of the galaxy. For this dependence, equa- 
tion ([4]) can be integrated to give 
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where R' = R/Rq, Rq = V m /u)Q and U' = 
U eS /(V^/2). For the case of Figures Q] and [2j 
V m m 555 km s -1 , lu q w 3.51 x 10" 15 s _1 , and 
Rq « 5.13 x 10 3 kpc. Of course, for V < V m , 
UcS = (l/2)wgi? 2 to a good approximation. Fig- 
ure [3] shows this function and its derivative which 
is the negative of the effective force, in contrast 
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Fig. 1. — The top panel shows the motion along 
the direction of kick (in the plane of the galaxy) 
for a BH of mass 10 8 M Q for an initial velocity 
kick of 500 km s _1 (solid line), 350 km s _1 (dashed 
line), and 150 km s _1 (dash-dot line). The bottom 
panel is for an initial kick velocity 200 km s _1 . 
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Fig. 2. — Amplitudes and frequencies of the cen- 
tral BH motion resulting from velocity kicks in the 
plane of the galaxy for a 10 8 M Q BH. 




Fig. 3. — Effective potential for the BH as it os- 
cillates about the galaxy center of mass after re- 
ceiving an impulsive kick. The effective restoring 
force is oc —dU' eS /dR'. The maximum of the force 
is « 0.36 at a distance R' s=s 0.56. 



to the purely logarithmic potential of the halo 
isothermal sphere. 

We have taken an example frame from the 
300 km s _1 model run, inclined it to the line of 
sight, and convolved the gas component with a 
Gaussian beam. The result is Figure [4] which de- 
picts the morphology and kinematics of the radio 
gas component of our galaxy model. Although 
the central BH is in motion, observable effects on 
the host galaxy are minimal. The simulated HI 
gas is distributed axisymmetrically, for instance. 
Moreover, as can be seen from the figure, lines of 
constant radial velocity are straight and parallel 
through the dynamical center of the galaxy, in- 
dicating pure rotation with little or no streaming 
motion, even though the central BH is clearly dis- 
placed from the dynamical center of the galaxy. 
This is the behavior we observed throughout this 
set of runs. 

A series of runs was also conducted with the 
velocity kick normal to the plane of the galaxy. 
Qualitative results were similar to the case where 
the kick is in the plane of the galaxy, with the BH 
oscillating roughly sinusoidally through the center 
of mass. The amplitude -initial kick relation re- 
mains linear for a slightly smaller range of kicks, 
however (Figure [5]) . For velocity kicks less than 
350 km s _1 , the amplitude-initial kick relation is: 
A(V )/kpc = (8.06 ± 0.51) x lO^Fo/kms" 1 + 
(0.126 ± 0.067), making the linear part of the 
amplitude-kick relation independent of direction 
to within 2a. Again, escape from the galaxy oc- 
curs for an initial kick 650 km s _1 . 

As in the coplanar case, little disturbance of 
the morphology and kinematics of the stellar or 
gas disk is observed; only a displacement of the 
BH from the center of the galaxy disk is ob- 
served. However, when viewed nearly face-on, 
some streaming motions are observed in the gas 
kinematics, as illustrated in Figure [6l 

3.2. Black Hole Mass 10 9 M Q 

When receiving an impulsive initial velocity 
"kick," a large 10 9 M Q black hole exhibits quali- 
tatively similar behavior to the 10 8 M Q black hole, 
as shown in Figure [7] In this case, however, it 
is clear that the damping time of the oscillation 
is significantly shorter. In the 10 8 M Q case, the 
damping scale time was on the order of 1.5 Gyr; 




Fig. 4. — Morphology and kinematics of the 
galaxy model at t = 0.12 Gyr after the central 
BH of mass 10 s M has been given a velocity kick 
of 300 km s _1 in the plane of the galaxy. In the 
top panel, the stellar disk is illustrated in greyscale 
and the logarithmic contour lines indicate density 
of HI gas, after convolution with a simulated Gaus- 
sian beam (lower left). In the lower panel, the 
convolved radial velocity of the gas is indicated by 
the contour lines, labeled in km s . The inclina- 
tion to the line of sight is 20°. In both diagrams, 
the BH is indicated by the black dot slightly to the 
right of center. Major tick marks delineate 25 kpc. 
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Fig. 5. — Amplitudes and frequencies of the BH 
motion resulting from velocity kicks perpendicu- 
lar to the plane of the galaxy for a BH of mass 
10 8 M n . 




Fig. 6. — Morphology and Kinematics of the 
galaxy model at t = 0.70 Gyr after the central 
BH has been given a velocity kick of 300 km s _1 
perpendicular to the plane of the galaxy for a BH 
of mass 10 8 M©. In the top panel, the stellar 
disk is illustrated in greyscale and the logarithmic 
contour lines indicate density of HI gas, after con- 
volution with a simulated Gaussian beam (lower 
left). In the bottom panel, the convolved radial ve- 
locity of the gas is indicated by the contour lines 
labeled in km s _1 . The inclination to the line of 
sight is 10°. In both panels the BH is indicated 
by the black dot in the center. Major tick marks 
delineate 25 kpc. 



however, for the 10 9 Mq case, damping times are 
of the order of 0.8 Gyr. 

During the period of significant motion of the 
black hole, significant dynamical disturbance can 
be detected in the HI velocity maps. After a 
damping scale time, however, the dynamics re- 
cover to a normal rotation curve. Morphology out- 
side the disk scale radius remains generally undis- 
turbed, although minor observable asymmetries 
occasionally appeared in the central gas distribu- 
tion. Despite the normal disk and gas morphology, 
however, as in the previous case, the central black 
hole could be found offset from the center of the 
disk by several kpc. 

Typical induced kinematic asymmetries are 
presented in Figure [8] In these plots, gas par- 
ticle velocities were convolved and combined in a 
simulated Gaussian telescope beam and lines of 
constant radial velocity drawn. If the galaxy were 
a symmetric rotator, then constant radial velocity 
lines should be straight and parallel as they pass 
through the dynamical center. However, as can 
be seen in the figure, these kinematics have been 
disturbed by the motion of the BH, which has re- 
sulted in a dynamical friction interaction with the 
galaxy. The central regions of the galaxy, within 
the disk scale radius of 3.5 kpc react quickly to 
the transferred momentum of the BH, and tend 
to move with it in its direction of motion. Beyond 
that radius, particles are shielded from the gravity 
of the central BH and react sluggishly. As a re- 
sult, the dynamical center of the central galaxy is 
displaced from the overall dynamical center in the 
direction of the initial impulse given to the BH. 
In the transition region between the two regimes, 
streaming motions are evident, with magnitudes 
on the order of 30 km s _1 . Finally, because of the 
BH's resulting orbit around the new dynamical 
center of the distribution, at most times its actual 
location was offset from the dynamical center as 
in the second frame of Figure [5] 

4. BH Acceleration Due to a One-Sided 
Jet 

In order to simulate the influence of a one-sided 
jet, we did a series of runs where the massive BH 
was placed in the center of the galaxy with zero 
initial velocity relative to the galaxy, but with 
an applied constant acceleration for a time of the 



order of the Salpeter time. The available accre- 
tion power from a BH is L a = eM a c 2 where M a 
is the mass accretion rate and the efficiency is 
assumed to be e = 0.1. Usually, the accretion 
rate is measured in units of the rate which would 
give the Eddington luminosity for the BH, that is 
M E dd = L Edd /(ec 2 ) » 2.2(AWlO 8 M Q )M yr-\ 
where L Edd » 1.26 x 10 46 (A/ 6ft /10 8 M Q ) erg/s and 
Mbh is the BH mass. Thus, the accretion lu- 
minosity can be written as L a — mL Edd , where 
rh = M a /M Edd . The e xponentia t ion ti me of the 
BH mass growth, the Salpeter ( 1964h time, is 
T s ~ 4.5 x lO 7 ^" 1 yr. 

It is widely thought that the jets of active 
galaxies are due to the presence of a strong, or- 
dered magnetic field (~ 10 3 ~ 4 G) threading the 
BH accretion disk. For the commonly considered 
case of a field with dipole symmetry the magnet- 
ically driven jets a re of equal stren gth on the two 



sides of the disk (Lovelace 1 119761 ). However, as 



IWang et al. I (|l992t ) pointed out, an ordered field 
in the disk consisting of dipole and quadrupole 
components can give magnetically driven jets of 
unequal strength on the two sides of the disk. In 
fact, one side of the disk may have a dominant jet. 
The timescale for the jet to be on one side r in this 
model is of the order of the time for plasma and 
magnetic field to move through the disk. This time 
may be shorter than Tg- Nevertheless, we empha- 
size in our simulations the influence of a one-sided 
jet by maintaining its force on the BH for a time of 
the order of Tg after which the force is zero. We 
assume that the BH carries along with it bound 
gas and stars sufficient to supply M a during this 
time. 

For the case of a one-sided jet, we can write 
the jet luminosity as i jot = fjetL a , where / jct < 1 
is the fraction of the accretion luminosity going 
into the jet. Alternatively, for asymmetrical, op- 
positely directed jets we take Lj Ct to be the differ- 
ence in the luminosities of the two jets. Assuming 
that the jet outflow is highly relativistic, the force 
on the BH/disk system is therefore 



F bh = -L jot /c« 4.2xl0 35 /jc t m ( 



V 1O 8 M 



dynes , 
(6) 



and the BH acceleration is 



a b h = -ji- 1 ~ 2.1 x 10 Vjetm -5- , (7) 

Mbh 
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(jWang et al. II1992T ) . The momentum imparted to 
the BH during a time T is I = L dt M bh a^h = 
TMbhO-bh for constant one-sided jet. For a jet 
which changes sides randomly on a timescale t <C 
T, the imparted momentum is much smaller, / ~ 
(rT) 1 / 2 : M bh a bh . The sm all values of the BH dis- 
placements estimated by IWang et al. I (|l992h re- 
sulted from an invalid assumption about the rigid- 
ity of the central potential of the galaxy. 

As a reference value for our simulations, we 
take ai,h = 2 x 10~ 8 cm s _1 which corresponds 
to fjetih — 0.01. For this value of abh, notice 
that dt a bh = Ta bh w 1260(T/2 x 10 8 yr) 
km s . We confirmed the relevance of this value 
by observing that the hole remains bound to the 
galaxy for applied accelerations up to a b h — 4.5 x 
10~ 8 cm s~ 2 , and escapes the galaxy for larger 
a b h- A hun dred times large r value of abh was con- 
sidered by iTsvgan (2007) to apply over a time 
scale Ts and this will strongly eject the BH from 
the galaxy. 

4.1. Black Hole Mass 10 8 M© 

When an external acceleration was applied in 
the plane of the galactic disk, the typical behavior 
was for the 10 8 M Q BH to find an equilibrium posi- 
tion off-center of the morphological and dynamical 
center of the galaxy and to oscillate about that po- 
sition as long as the acceleration was applied. Fig- 
ure[9]illustrates the BH displacement as a function 
of time for two typical examples of this experi- 
ment, in which the 1O 8 M black hole was given 
an applied acceleration of 2.0 x 10~ 8 cm s" 1 and 
4.0 x 10 -8 cm s _1 in the plane of the galaxy. The 
center of mass of the galaxy is also depicted for 
the first experiment. This accelerates as well, as 
dynamical friction transfers momentum between 
BH and the galaxy. In the second run, the BH es- 
capes the galaxy entirely within 0.1 Gyr. In each 
run where the BH remained bound, we confirmed 
that the acceleration of the center of mass of the 
galaxy was equal to abh(Mbh / M to t) , where abh was 
the acceleration imparted to the BH, and also that 
when abh was turned off after 200 Myr, that the 
final velocity of the center of mass agreed with the 
total impulse I = J^ s dtMbhO-bh given, within er- 
rors which were taken to be the observed fluctua- 
tion in the center of mass velocity and acceleration 
of a galaxy with no kick (approximately 10%). 



For BHs of mass 10 8 Af Q , slight to moderate 
asymmetries were induced in the dynamics of the 
host galaxy during the "on" phase, although noth- 
ing detectable was induced in the morphology, as 
illustrated in Figure [TO] When viewed with the 
galaxy inclined to the line of sight, with the axis 
inclination aligned with the applied acceleration, 
deformations of the contours of constant radial 
velocity of order 40 km s^ 1 are clearly observed 
in the kinematics of the galaxy gas. Note that 
no grand-design spiral arms have formed in the 
stellar or gaseous components, ruling out density 
waves as a possible explanation for these. These 
streaming motions are also present when the line 
of sight is at other orientations to the line of sight, 
but are less prominent. These disk asymmetries 
disappeared when the acceleration was turned off, 
but the BH continued to oscillate about the center 
of mass with an amplitude that depended slightly 
on the initial acceleration, but was typically of or- 
der 1 — 2 kpc. 

In the case where the BH is given an accelera- 
tion perpendicular to the galaxy plane, a similar 
behavior is observed, as illustrated in Figure [11] 
The BH oscillates around an equilibrium point, 
slowly dragging the galaxy along with it via dy- 
namical friction. In this case, however, the ampli- 
tude of the motion is smaller and the frequency 
higher. The acceleration required for escape in 
this case is 3. Ox 10~ 8 cm s -2 , as compared with 4.5 
for the case in the plane of the galaxy. This case 
produces no detectably significant asymmetries in 
the morphology or kinematics of the galaxy disk. 

4.2. Black Hole Mass 1O 9 M 

When given an acceleration in the plane of the 
galactic disk, the large 10 9 M Q black hole exhib- 
ited behavior similar to that of the smaller hole. In 
this case, as illustrated in Figure IT2l the larger im- 
pulse "drags" the galaxy a proportionally greater 
distance. When the acceleration is turned off af- 
ter 0.2 Gyr, the velocity imparted to the galaxy is 
200 m s -1 . 

As the black hole was accelerated, it became 
displaced from the morphological and dynamical 
center of the galaxy, as illustrated in Figure [131 
These induced asymmetries in the galaxy disk are 
more pronounced in this case than in the 1O 8 M0 
case, in particular tending to drag the central re- 
gions of the galaxy with the movement of the 



8 



black hole via dynamical friction. In particular, 
while the BH in the first frame of Figure [13] ap- 
pears displaced from the morphological center of 
the galaxy, the second frame of the figure illus- 
trates how the dynamical center of the galaxy has 
been displaced towards the direction of accelera- 
tion of the black hole. The result is a character- 
istic "tongue-" shaped extension of the velocity 
contours on the side of the galaxy opposite the 
acceleration and flattened contours on the side of 
the galaxy in the direction of the acceleration. In 
test runs with a much larger, possibly unphysi- 
cal 1O 1O M BH, these features became associated 
with a central core region of the disk that moved 
along with the black hole, being more coupled to 
it than to the galaxy, and causing strong mor- 
phological asymmetries in the disk as the bulge 
was dragged off-center. We did not observe such 
strong asymmetric features in the morphology of 
the 10 9 Af Q galaxy, but do observe the "tongue" 
features that indicate that the BH is beginning to 
dominate the dynamics of the inner 2 — 5 kpc of 
the disk. 

A similar behavior is observed when the acceler- 
ation is perpendicular to the plane of the galaxy. 
Again, streaming motions and dynamical asym- 
metries are observed for the first 0.5 Gyr of ac- 
celeration, as shown in Figure [14] Following that 
period, only slight streaming motions remain as 
part of the stabilized disk kinematics. 

5. Discussion 

We have analyzed the motions of massive black 
holes in the centers of galaxies which have received 
an impulsive kick (owing to the merging of two two 
BHs) or received a sustained acceleration due to a 
one sided jet. We have studied the BH displace- 
ment from the galaxy center and the influence this 
has on the galaxy morphology and kinematics. An 
impulsive "kick" exceeding 600 km s _1 will ejects 
the BH from the galaxy. Accelerations of the order 
of 4 x 10~ 8 cm s -2 or larger for times of the order 
of 2 x 10 8 yr eject the BH from the galaxy. This 
acceleration is smaller than th e Eddington limit 
on the one-sided jet luminosity ( Tsvgan 120071 ) by 
a factor of about 50. 

For smaller velocity kicks or smaller sustained 
accelerations, the BH is observed to oscillate 
through the galaxy center for times of order Gyr. 



The frequency of these oscillations decreases as 
the BH velocity kick increases or as the BH accel- 
eration increases. In such galaxies, we can expect 
to observe central BHs which are offset from the 
centers of their host galaxy disks or have large 
radial velocity differences from their hosts. Larger 
mass BHs undergoing these oscillations induce 
streaming motions in their galaxies due to dy- 
namical friction. This process may be a source of 
disk morphological and dynamical asymmetry in 
galaxies such as NGC 1637 or NGC 991 which ex- 
hibit nonaxisymmetries but are located in isolated 
fields far from potential sou rces of tidal interac- 



tions ( Kornreich et al. 19981 ) 



The asymmetries in the galaxy disk induced by 
the BH motion are of course more pronounced for 
more massive black holes. The BH tends to drag 
the central regions of the galaxy with it owing to 
dynamical friction. The dynamical center of the 
galaxy is seen some cases to be displaced towards 
the direction of the BH acceleration with a charac- 
teristic "tongue-" shaped extension of the velocity 
contours on the side of the galaxy opposite the ac- 
celeration. 

As a check on the results, additional model runs 
with fewer (10 5 ) or more (10 7 ) particles yielded 
substantially similar behavior for the BH. Asym- 
metries induced in the galaxy disk remained qual- 
itatively similar for each run, although the de- 
tails of the locations and intensities of the induced 
asymmetries varied somewhat. 

It is important to extend these studies to ellip- 
tical galaxies in the future. We expect that the 
results concerning conditions of BH escape from 
the galaxy will apply equally well there, as es- 
cape times from the disk are very short, of order 
0.1 Gyr, and after that the gravitational potential 
is essentially that of the spheroidal halo. 

We thank Martha Haynes and Larry Kidder for 
discussions. This work has made use of the compu- 
tational facilities of the National Astronomy and 
Ionosphere Center, which is operated by Cornell 
University under a cooperative agreement with the 
National Science Foundation. 
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Fig. 7. — Motion along the direction of kick, in 
the plane of the galaxy, for the Black hole of mass 
10 9 M© inside a model disk galaxy for an initial 
velocity kick of 400 km s -1 (solid line), 300 km s _1 
(dashed line), and 100 km s _1 (dotted line). 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 
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Fig. 8. — Kinematic asymmetries in the gas com- 
ponent of a galaxy with a 1O 9 M BH which has 
received an impulsive velocity kick. In the top 
panel, the kick velocity is 200 km s _1 to the right 
in the plane of the galaxy. The galaxy is inclined 
to the line of sight by 30°. The bottom panel is 
for the same parameters but for an initial velocity 
kick of 400 km s _1 . Contours of constant radial 
velocity are labeled in km s _1 . Major tick marks 
delineate 25 kpc. 
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Fig. 9. — Displacement vs. Time along the di- 
rection of applied acceleration in the plane of the 
galaxy for a 10 s M Q black hole with simulated 
jet acceleration of 2.0 x 10~ 8 cm s~ 2 for a time 
< t < 0.2 Gyr. The dotted line shows the 
center of mass of the galaxy. The bottom panel 
shows the same quantities but for an acceleration 
4.0 x 10~ 8 cm s~ 2 . 
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Fig. 10. — Galaxy model gas density contours and 
velocity map, convolved with a simulated Gaus- 
sian telescope beam, for the 1O 8 M0 BH with an 
applied acceleration of 2.0 x 10~ 8 cm s~ 2 at time 
t = 0.2 Gyr. The inclination to line of sight is 30° 
about the axis of acceleration. Major tick marks 
delineate 25 kpc. 



2.5 r 




-0.5 L_, , , I , , , I , , , I , , , I , , ^_ 

0.0 0.2 0.4 0.6 0.B 1.0 

Time (Gyr) 

Fig. 11. — The BH displacement as a function of 
time along the direction of the applied accelera- 
tion perpendicular to the plane of the galaxy for a 
10 8 M© black hole. The simulated jet acceleration 
is 2.0 x 10 ~ 8 cm s _1 . The dotted line shows the 
center of mass of the galaxy. 
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Fig. 12. — Displacement vs. Time along the di- 
rection of applied acceleration for the black hole 
of mass 10 9 M© with simulated jet acceleration of 
2.0 x 10~ 8 cm s~ 2 . Acceleration in the plane of 
the galaxy is represented by the solid line; per- 
pendicular to the plane by the dotted line. In the 
bottom panel, the acceleration has been shut off 
at 0.2 Gyr. 



Fig. 13. — Galaxy model gas density contours 
(top panel) and velocity map (bottom panel) , con- 
volved with a simulated Gaussian telescope beam, 
for the case of a 10 9 M Q black hole with an applied 
acceleration of 2.0 x 10~ 8 cm s~ 2 in the plane of 
the galaxy, at time t — 0.2 Gyr. The inclination to 
line of sight is 30° about the axis of acceleration. 
Major tick marks delineate 25 kpc. 
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Fig. 14. — Gas density contours (top panel) and 
velocity map (bottom panel), convolved with a 
simulated Gaussian telescope beam, for the case 
of a 10 9 M Q black hole with an applied accelera- 
tion of 2.0 x 10~ 8 cm s _1 , perpendicular to the 
plane of the galaxy at time t = 0.2 Gyr. The in- 
clination to the line of sight is 30°. Major tick 
marks delineate 25 kpc. 
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